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Abstract 

The  estimation  of  individual  instabilities  of  N  clocks,  when  only  differences  of  phase-time 
clocks  readings  are  available,  can  be  carried  out  even  supposing  cross-correlation  between 
clocks.  Based  on  the  original  analysis  of  this  problem,  developed  by  P.  Tavella  and  A. 

Premoli,  and  the  Kuhn-Tucker  theorem  (1951),  a  new  way  to  solve  the  constrained 
minimization  problem  will  be  introduced.  We  propose  our  particular  Kuhn-Tucker 
equations:  The  already  well-known  constraint  function  and  two  new  candidates  objective 
functions.  Advantages  and  inconveniences  in  the  objective  function  election  will  be 
considered.  Finally,  our  proposal  will  consist  of  solving  two  constrained  minimization 
problems:  the  first  of  them  to  obtain  an  intermediate  solution  to  be  used  as  initial  condition 
in  the  second  minimization  problem.  Examples  with  some  experimental  data  from  the  Real 
Observatorio  de  la  Armada  en  San  Fernando  illustrate  the  capabilities  of  this  proposal. 

INTRODUCTION 

It  is  well  known  that  an  useful  mathematical  tool  to  characterize  the  stability  of  any  device  is  an 
estimation  of  the  variance  by  using  the  available  measurements.  When  the  stability  of  several  clocks 
is  evaluated,  the  absolute  variance  cannot  be  directly  estimated  because  the  measurement  data 
available  are  time  deviations  between  pairs  of  real  clocks. 

In  past  years,  many  methods  were  introduced  to  solve  this  problem:  from  the  classical  and  popular 
“three-cornered-hat”  method,  to  the  subsequent  generalizations  to  N  clocks  [1],  but  all  the  works 
assume  the  hypothesis  of  independence  between  clocks  (and  therefore,  the  hypothesis  of  null  cross¬ 
correlation  variables).  The  problem  so  resolved  occasionally  produces  negative  estimates  of  variance; 
this  fact  may  be  attributed  to  several  causes,  the  more  notable  one  to  accept  incorrelation  between 
clocks. 

In  a  recent  work,  A.  Premoli  and  P.  Tavella  [2]  expound  lifting  the  assumption  of  uncorrelation,  and 
proposed  a  revisited  version  of  the  three-cornered-hat  method,  consistent  and  formally  equivalent  to 
the  classical  one  when  clocks  are  uncorrelated,  but  that  no  longer  produces  negative  estimates  of 
variance.  Later,  the  same  authors  establish  the  theoretical  basis  to  undertake  the  problem  generalized 
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to  N  clocks,  and  demonstrate  that  the  arbitrariness  in  estimating  the  solution  is  reduced  when  the 
number  of  compared  clock  increases  [3]. 

At  30th  PTTI,  F.  Torcaso,  C.R.  Ekstrom,  E.A.  Burt,  and  D.N.  Matsakis  proposed  a  numerical  method 
for  the  practical  resolution  of  the  problem  generalized  to  N  clocks  [4], 

The  goal  of  this  work  is  to  propose  an  alternative  method,  based  on  the  use  of  a  new  objective 
function,  with  a  better  physical  meaning.  The  technique  will  be  applied  to  simulated  and  real  series  of 
data.  The  obtained  results  will  be  finally  discussed. 

NOTATION  AND  DEFINITIONS 

The  absolute  estimated  variance  can  only  be  calculated  from  time  deviations  between  a  physical 
clock  and  an  ideal  one,  but  a  good  approach  can  be  carried  out  from  the  existing  relations  between 
absolute  variances  and  covariances  (jointly  called  (co)variances)  and  their  corresponding  relatives 
between  individual  clocks. 

Let  X '  =  {x [  :  k  >  l}  for  i  =  N  be  the  time  differences  of  the  t*  clock  referenced  to  an  ideal 
one,  measured  at  intervals  of  T0  seconds.  Denoting  the  absolute  fractional  frequency  of  the  j*  clock, 
averaged  over  a  time  Z  =  mt0  starting  at  time  tk  ,  by: 


n'W= 


a) 


The  Allan  variance,  or  two-sample  variance,  can  be  expressed  as: 

(2) 

where  Y, ,  Y2  denote  pairs  of  adjacent  absolute  fractional  frequencies  and  (  )  denote  mathematical 

expectation  or  infinite  time  average.  We  suppose  here,  as  usual,  that  the  process  of  averaged  fractional 
frequencies  is  stationary  and  ergodic,  making  the  expression  (2)  well-defined. 

Let  now  x 1  =  \x‘k  :  k  >  l}  for  i  =  1, . . . ,  N  —  1  represent  the  time  differences  of  the  Ith  clock  referred  to 

the  A*  one  taken  as  a  reference;  then  we  can  express  x'  =  X'  —  X‘v .  The  corresponding  fractional 
frequencies  are: 


Y*(T)  = 


x 


k+m 


z 


(3) 
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The  absolute  and  relative  fractional  freauencies  are  related  bv  the  following  expression- 

A  •>  W  jC - 

y[  (t)  =  Yk  -  YkN  ,k>l,  and  the  Allan  variance  of  i*  clock  referenced  to  A*  clock  can  be  expressed 

^■W=5((yi-y;)!)  (4) 

by: 

Both  su  (t)  and  r  (t)  are  theoretical  quantities  impossible  to  obtain  experimentally.  They  could  only 
be  estimated  from  a  finite  number  of  samples  yk  (r)  .  The  estimate  of  Allan  variance  can  be  obtained 
by  the  expression: 


1 


/-s  \ 


z(M  -zmjm  r o 


M -2m 

X  ( xk+2m-2xk+m+xk j 


(5) 


In  this  paper,  the  Allan  covariance  of  clocks  Ith  and  _/*  referenced  to  the  A*  clock  are  needed, 
moreover,  we  will  attempt  to  estimate  the  dereferenced  Allan  covariances  between  pair  of  clocks.  We 
define  both  statistical  parameters  as  they  were  introduced  in  [4]: 


y,M=i/(v;-y;  (6) 


We  introduce  again  the  estimates  of  both  quantities  as: 


»  IJ\. 
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It  is  obvious  that  a  covariance  is  reduced  to  a  variance  making  /=/,  so  it  is  easy  to  derive  the 
expression  for  f  from  A  . 


In  the  subsequent  paragraphs,  we  will  assume  T  as  the  time  interval  in  which  we  want  to  characterize 
the  stability,  so  we  will  omit  the  symbol  T  assuming  that  time  interval. 

According  to  the  previous  expressions,  the  following  relations  can  be  deduced: 

c  —  *■  _L_  r  _  *-  y*  YQ\ 

0  ij  '  ij  '  1 NN  riN  ‘  jN  V*-V 


and 
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These  quantities  can  be  jointly  denoted  as  matrices.  The  N  by  N  matrices  R  and  R  are  related  to  the 
dereferenced  (co) variances;  on  the  other  hand,  the  (N- 1)  by  (N- 1 )  matrices  S  and  $  are  related  to  the 
(co)variances  referenced  to  a  Nth  clock. 

THE  PREVIOUS  APROACHES 

Before  introducing  the  new  proposal  of  this  work,  we  need  to  make  a  review  on  the  previous 
approaches  based  on  similar  arguments,  that  is:  the  Tavella-Premoli  approach  [2],  [3]  and  the  F. 
Torcaso  et  al.  approach  [4], 

In  a  first  approach,  Tavella  and  Premoli  propose  rejecting  the  hypothesis  of  uncorrelated  clocks, 
leading  to  an  underdetermined  linear  system.  They  introduce  the  estimates  of  covariance  Sl}  and  a 

suitable  optimization  criterion  to  calculate  the  complete  covariance  matrix  R  .  Then,  they  analyze  in 
detail  the  case  of  three  clocks,  obtaining  finally  a  quasi-analytical  solution. 

According  to  their  definitions,  both  $  and  R  are  symmetric  matrices,  so  the  number  of  unknown 
quantities  is  reduced  to  iV(iV  +  l)/2,  while  the  number  of  equations  is  {N  —  l)N/2.  For  the  three- 
clock  case,  the  six  unknown  quantities  are  the  (co)variances  rn ,  r22 ,  r33 ,  rl2 ,r]3>  and  r23  and  the 

three  equations  are: 

$11  =  *11  ^33  —  2^13  (10) 

^12  =  **12  +  *33  _  **13  _  *23 
$22  ~  $22  +  *33  ~  “  $22 

Nevertheless,  they  outlined  an  important  constraint  which  bounds  the  solution  space  and  guarantees  a 
significant  result:  the  positive  definiteness  of  any  covariance  matrix.  They  expressed  this  constraint  as 
a  function  H  (f13 ,  f23 ,  f33 ),  so  that  the  matrix  R  would  be  positive-definite,  provided  that  the  function 

H{$\3  >$22  >  *33  )  >  0  : 

H{$\3>$23>$33)=  $33  ~($12  ~  $33’$22  ~  $3l)S  (*1 3  —  *33  ’  *23  ~  *33  )  (1  1) 


This  function,  when  equalled  to  zero,  represents  the  bound  of  the  solution  domain,  and  geometrically 
it  represents  an  elliptic  paraboloid.  Afterwards,  the  authors  suggested  the  determination  of  the  values 
r13 ,  r23,and  r33  by  means  of  the  resolution  of  a  minimization  problem,  and  they  defined  an  objective 

function  [G($n ,  $23 ,  f33  )]2that  is  directly  proportional  to  the  quadratic  sum  of  covariances  and 
inversely  to  H: 


F{$\3>  $23  ’$ 33  ) —  K 


[<j(fi3  ,  $23  ,  f33  )] 

H  (f,3  ,  $22  ,  $33  ) 


(12) 
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Therefore,  they  minimized  the  “global  correlation”  among  clocks  maintaining  the  positive  definiteness 
of  R  leading  to  the  values  r13,  r23,and  r33.  The  quantities  that  minimize  the  objective  function  are 

then  substituted  in  (10)  concluding  the  calculation  of  R  . 

Because  of  the  convex  character  of  the  objective  function,  the  solution  to  the  minimization  problem  is 
unique.  This  solution  is  reduced  to  the  classical  three-cornered- hat  solution  provided  that 
5,,  >Sl2,S22  >512  and 5I2  >0. 


In  a  second  paper,  Tavella  and  Premoli  look  into  the  generic  problem  for  N  clocks.  They  demonstrate 
that  the  matrix  R  is  positive  definite  if  and  only  if  R  >0,  deriving  a  compact  expression  to  define 

the  constraint: 


»('i. 


~  VNN  VI N  ^NN  ’ '  '  '  »  ?(N-l)N  ^NN  V  IV  ^AW  >•"’  ^(V-l)V 

whose  bound  can  be  determined  by  making  H(pin  fNN )  =  0 .  This 
geometrically  an  elliptic  hyper-paraboloid. 


(13) 

bound  represents 


Another  interesting  conclusion  formulated  is  that  the  domain  delimited  by  the  comparison  of  N  clocks 
has  a  higher  dimension  due  to  the  additional  variable  ,  but  that  for  any  value  of  this  parameter, 

it  can  never  be  larger  than  the  domain  obtained  from  the  comparison  of  AM  clocks. 

The  reduction  in  the  amount  of  arbitrariness  in  the  determination  of  the  covariance  matrix  R  led  F. 
Torcaso  et  al.  [4],  to  think  about  the  generalization  of  the  Tavella-Premoli  scheme. 


These  authors  suggest  to  modify  the  objective  function  (12).  The  new  objective  function,  denoted  as 
the  Tavella-Premoli  function,  is: 


H  (flN  ,...,PNN) 


(14) 


where  H  is  given  by  (13).  The  presence  of  H  squared  in  the  denominator  keeps  the  minimization 
problem  scale-invariant  and  facilitates  the  numeric  resolution.  The  initial  proposed  conditions  are 
selected  in  order  to  assure  that  the  initial  data  lie  between  the  constraint: 

f“'=0for  i<N  (15) 

C=4t.  *'=(1 . l)S-‘  (I - l)r 

2  s 
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The  minimization  problem  has  a  unique  solution  due  to  the  considerations  of  convexity  already 
mentioned.  The  N  variables  ftN  =  i  =  1,..., N ,  can  be  determined  this  way.  The  remaining 
quantities  are  determined  by  resolving  the  system  of  (N  —1)N/2  equations  stated  in  (9). 

THE  CONSTRAINED  MINIMIZATION  PROBLEM 


The  reduction  of  arbitrariness  on  the  solution  domain  due  to  the  increasing  of  the  number  of  clocks  is 
a  justified  reason  in  order  to  look  for  other  valid  alternatives  that  generalize  the  expounded  problem 
by  Tavella  and  Premoli. 

The  goal  of  this  new  proposal  is  to  reach  a  solution  in  close  agreement  with  the  physical  reality.  This 
can  be  attained  avoiding  the  distortion  of  the  objective  function  with  the  constraint  condition.  This  is 
the  reason  we  have  used  the  Kuhn-Tucker  (AT)  theorem,  published  in  1951.  This  theorem  states  the 
following: 

“ Given  a  problem  of  the  type: 

min  /(x)  (16) 

g,(x)<0,  i  =  l,...,m 


with  a  local  minimum  in  Xmin  .  If  the  gradient  vectors  of  the  saturated  restrictions  in  xmin  are  linearly 
independent  over  Xmin ,  then  each  saturated  restriction  has  a  positive  number  \  >  0  (known  as  KT 
multiplier)  associated,  that  verifies: 

V/(x.,J+tA(-Vg,(x.J  =  0  (17) 

1=1 

so  that,  in  the  previous  expression,  all  the  null  Multipliers  appearing  in  the  restrictions  are  associated 
to  the  not  saturated  restrictions.  In  order  to  get  this,  the  following  conditions  are  imposed: 


V<?,(O  =  0>  i  =  l,...,m. 


(18) 


A  restriction  is  said  to  be  saturated  when  the  optimal  solution  is  verified  as  equality,  otherwise  the 
restriction  is  said  not  saturated. 

The  KT  equations  are  necessary  conditions  for  optimality  for  a  constrained  minimization  problem.  If 
the  problem  is  a  so-called  convex  programming  problem,  that  is,  / (x)  and  g;  (x)  i  =  1, . . . , m.  are 

convex  functions,  then  the  KT  equations  are  both  necessary  and  sufficient  conditions  for  a  global 
solution  point.  If  the  objective  function  is  also  strictly  convex,  in  the  case  the  minimum  exists  that  is 
unique. 

The  KT  theorem  suggests  the  employment  of  a  constraint  function  derived  from  the  (13)  condition. 
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modified  appropriately  so  that  make  possible  the  application  of  the  KT  theorem: 

0  (19) 

K 

The  strict  inequality  condition  that  comes  in  (19)  could  seem  an  important  limitation,  but  it  is  not  so 
from  a  practical  point  of  view.  On  the  other  hand,  K  denotes  K  -  has  been  introduced  for  the 

sake  of  adimensionality  and  to  facilitate  the  numeric  resolution.  This  constraint  function  represents 
geometrically  the  interior  of  an  elliptic  hyper-paraboloid,  as  was  seen  previously,  being  therefore  a 
convex  function. 


Two  objective  functions  could  be  candidates  in  the  minimization  problem,  each  one  with  their 
advantages  and  disadvantages.  The  first  one  seems  to  be  obvious,  knowing  the  antecedents  of  the 
problem: 


i^\N  »•  ^AW  ) 


y  k2 

£—ii<  j  u 


(20) 


The  second  one  has  greater  physical  meaning,  although  it  presents  certain  inconveniences  when  it  is 
used  in  a  minimization  problem: 


^2  i^lN 


'  NN 


)-L 


i<  j  P 
ii 


(21) 


The  convexity  of  the  function  done  by  (20)  could  be  easily  demonstrated  through  the  Hessian  matrix 
because  the  second  derivatives  are  constant.  A  function  that  admits  derivation  partially  twice  is 
convex  if  and  only  if  its  Hessian  matrix  is  positive-semidefinite  for  any  value  in  the  domain.  If  the 
Hessian  matrix  is  positive-definite,  the  function  will  be  strictly  convex.  The  obtained  Hessian  matrix 
eigenvalues  are  the  following: 

Xi=2{N-2\  i  =  \,...,N -2  (22) 

N2  +N -4-yjN4  +2N3  -15N2  +16N 
Xn~1  ~  2 

N2  r  N  -4-\  -  N4  +2N3 -15N2  +16N 


For  N  >  3 ,  we  conclude  that  Fx  function  is  strictly  convex. 

Regarding  to  the  second  alternative:  F2,  it  could  not  be  sure  that  an  unique  solution  exists.  We 

illustrated  an  example  in  Figure  1  in  which  this  fact  is  stated.  In  that  Figure  are  shown  slices  of  the 
elliptic  paraboloid,  that  limits  the  solution  domain  for  S  =  (l. 09, 1.18;  1.18, 11.35).  S  has  been 

calculated  over  three  simulated  series  of  gaussian,  white  frequency  noise  data,  with  Allan  variances 
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and  correlated  coefficients  predetermined:  tu  =  2,t22  =  16,f33  =  1  y  j  =  0.68  Vi*  y  .  Figures  la,  lb 

and  lc  shown  the  planes  f33  =  (2-(l  l)$~‘  (l  l)r )  '  =0.54,  f33  =  0.75,  and  f33=1.00  including  the 
initial  point  (0,0,0.54)  and  the  prospective  solution  (0.97,2.73,1.00) .  It  can  be  seen  in  each  slice  the 
regions  where  F2  is  positive-definite  (+)  or  it  is  not  (-).  It  seems  to  be  obvious  that  the  function  (21) 

does  not  lead  to  an  unique  solution.  In  any  case,  this  function  behaves  in  a  soft  way,  without  wide 
variations. 

The  last  consideration  is  the  main  reason  why  a  good  election  of  the.  initial  conditions  is  considered 
essential:  The  closer  the  initial  condition  to  the  final  solution  is,  the  higher  will  be  the  probability  to 
converge  to  it. 


The  final  suggestion  consists  of  the  following:  In  a  first  approximation,  equations  (20)  and  (19),  with 
the  initial  conditions  (15),  allow  calculation  of  a  unique  solution  (minimum  of  the  quadratic  mean 
of  covariances).  This  solution  is  used  as  initial  condition  over  N  minimization  problems,  each  one 
done  by  a  rotation  of  the  N  clocks,  taking  in  each  case  a  different  clock  as  reference  clock.  The  value 
reached  by  F2  for  the  final  solution  should  be  lower  than  the  initial  value,  so  the  N  possible  solutions 

are  better  than  the  solution  obtained  in  the  first  step  (from  a  global  minimum  of  the  quadratic  mean  of 
correlation  coefficients  point  of  view).  The  last  step  consist  of  selecting  some  of  the  N  solutions;  for 
this,  we  could  apply  some  validation  test  able  to  select  the  best  solution.  The  used  test  is  based  on  two 
criteria: 


1 .  Repetitiousness  of  each  certain  final  solution 

2.  Homogeneity  of  the  obtained  cross-correlation  coefficients. 


The  first  criterion  is  the  more  important,  its  application  being  enough.  When  the  first  criterion  fails, 
the  second  one,  that  selects  the  solution  that  makes  smaller  and  more  homogeneous  the  absolute 
values  of  the  correlation  coefficients,  is  applied: 


min 


\i  {N  -  \)N  ^ i<j 


L<>,!-p)2 


with  p  = 


2 

{N-\)N 


(23) 


The  robustness  of  the  algorithm  could  increase  considering  the  N  solutions  reached  when  initial 
conditions  (15)  are  applied  minimizing  F2 . 

If  there  are  any  suspicions  that  two  clocks  are  more  or  less  correlated  than  the  others,  F}  and  F2 
might  be  formulated  by  introducing  a  covariance  factor  for  each  covariance  term. 
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RESULTS 


To  illustrate  the  technique  described  in  the  previous  section,  simulated  series  of  data  with  a  pre-set 
covariance  matrix  R  were  used,  having  the  considered  cross-correlation  \pu  =0.10  Vi*  j  . 

We  will  pay  attention  over  the  individual  variances  specially  (diagonal  of  R ). 


When  three  clocks  were  considered,  the  results  obtained  were  better  than  those  obtained  with  the  two 
previous  approaches,  but  we  have  a  relative  error  that  comes  to  be  up  to  75%.  The  solution  improves 
drastically  when  four  or  more  clocks  are  considered.  This  can  be  seen  in  Table  1,  where  the  results 
obtained  by  means  of  this  method  and  by  the  previous  one  proposed  by  F.  Torcaso  et  al.,  are 
indicated,  after  working  the  covariance  matrix  (24)  calculated  from  (7). 


^2.78 

0.95 

2.10  ' 

$= 

0.95 

4.60 

2.58 

2.10 

2.58 

394.57 

fu 

hi 

True  value 

2.03 

4.03 

Estimated  by  F.  Torcaso.  et  al.  method. 

2.09 

3.22 

New  proposal:  Intermediate  solution. 

1.60 

2.95 

New  proposal:  Final  estimate. 

1.81 

3.64 

Table  1.  Comparison  among  true  and  estimated  Allan  variances. 


(24) 


hi 

*44 

397.44 

1.02 

391.43 

2.14 

391.77 

1.88 

392.45 

0.97 

Table  2  shows  the  relative  error  made  for  each  case. 


clock  1 

clock  2 

clock  3 

clock  4 

Estimated  by  F.  Torcaso  et  al.  method. 

0.03 

0.20 

0.02 

1.10 

New  proposal:  Intermediate  solution. 

0.21 

0.27 

0.01 

0.84 

New  proposal:  Final  estimate. 

0.11 

0.10 

0.01 

0.05 

Table  2.  Relative  errors  ( ). 

XTrue 


Tables  1  and  2  don’t  show  significant  differences  between  the  estimate  proposed  by  F.  Torcaso  et  al. 
and  our  intermediate  solution,  but  both  estimates  introduce  high  relative  errors:  the  lower  the  Allan 
variance  is  the  higher  is  its  relative  error.  These  inadequacies  don’t  seem  to  exist  in  our  final  solution. 
Nevertheless,  we  should  keep  in  mind  that  the  simulated  data  used  to  calculate  the  matrix  S  have 
equal  cross-correlation  values  (in  absolute  value),  which  could  explain  the  good  results  obtained. 
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35  718  14  896  16113  121223  141569  31422  35  583 


Referenced  to  UTC  scale  0. 14 

Estimated  by  F.  Torcaso.  et  al.  method.  1.00  : 

New  proposal:  Intermediate  solution.  1.93  1.79  35.75  100.83  6.62  7.37  1 

New  proposal:  Final  estimate.  0.20  1.58  30.25  93.34  8.65  6.28  0 


Table  3.  Comparison  of  estimated  Allan  variances  evaluated  according  to  several  alternatives. 


The  technique  has  also  been  applied  on  phase-time  deviation  data  obtained  from  atomic  clocks  at  the 
Real  Observatorio  de  la  Armada  en  San  Fernando  (ROA).  Data  correspond  to  seven  commercial 
cesium-beam  standards  over  a  period  of  two  years  (1996-1997):  standards  583  and  718  are  HP  5071 
high  performance  models;  the  896,  569  and  1569  are  HP  5061  Opt. -004;  the  1223  was  an  HP  5061 A 
and  the  113  was  an  OSCILLOQUARTZ  model  3200;  the  series  of  data  were  taken  at  5-day  intervals, 
with  the  purpose  of  knowing  the  phase-time  deviation  relative  to  UTC.  The  Allan  covariance  matrix 
was  evaluated  for  an  integration  time  of  20  days,  with  overlapping  samples;  results  were  compared 
using  the  Allan  covariance  matrix  related  to  the  (co)variances  referred  to  the  UTC  scale  (for  a  time 

_ i  —  on  j  _  _ :  j _ i  Ai.  _  A  Ai.  _  t  rm _ i  _  •_  _  _ •  ^  .  . .  _  .  _  _  i_  i  _  _  .  ...  .  r _ _  _  i  _  _  i . .  .  .i 

interval  i  =  z,uu.,  we  cunsiuereu  urai  uie  uil  scaie  is  quite  more  siaoie  man  any  01  our  ciocks,  anu 

thus  could  be  supposed  quasi-ideal).  The  selected  time  interval  is  considered  wide  enough  so  that 
correlation  exists  among  the  clocks.  The  obtained  results  are  shown  in  Table  3,  where  numerical  values 
have  been  scaled  by  1028: 

At  this  point,  we  would  like  to  emphasize  that  the  time  period  was  selected  in  this  way,  because  in  this 
epoch  old  and  new  clocks  cohabited. 

When  we  look  at  the  Table  3  thoroughly,  two  results  stand  out  from  the  rest: 

1.  No  prominent  difference  exists  between  the  procedure  described  by  F.  Torcaso  et  al.  [4],  and  the 
intermediate  solution  obtained  when  the  constraint  condition  is  separated  from  the  objective 

fnnrtiAn 

±.  UilVUUll. 

2.  Results  estimated  with  the  new  proposal  seem  to  be  more  in  agreement  with  those  “expected”: 
usually,  minimizing  the  quadratic  sum  of  covariances  produces  a  pessimistic  estimate  of  the 
stability  for  the  most  stable  clocks  when  there  is  strong  heterogeneity  among  thorny  this  doesn't 
happen  with  the  new  proposal. 

CONCLUSIONS 


Separating  the  constraint  condition  of  the  objective  function  doesn't  seem  to  introduce  important 
benefits,  because  the  solutions  reached  in  both  cases  are  very  similar.  Nevertheless,  this  new  estimate 
truly  makes  minimum  the  quadratic  sum  of  covariances. 


The  resolution  of  the  problem  in  two  phases  supposes  an  important  qualitative  improvement  in  the 
estimation  of  the  solution,  since  the  second  objective  function  used  has  greater  physical  meaning.  The 
existence  of  several  minima  for  this  objective  function  should  not  be  a  problem  when  applying  the 
procedure  that  has  been  described  above. 

The  minimization  algorithm,  as  well  as  the  validation  tests,  is  easy  to  implement,  which,  along  with 
next  to  the  goodness  of  the  obtained  results,  could  convert  it  in  a  very  useful  method  for  estimating 
the  frequency  stability. 
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